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Abstract 

The propagation of a virtual quark in a thermal medium is considered. The non-perturbative 
jet transport coefficient q is estimated in quark less SU(2) lattice gauge theory. The light like 
correlator which defines q, defined in the regime where the jet has small virtuality compared to 
its energy, is analytically related to a series of local operators in the deep Euclidean region, where 
the jet's virtuality is of the same order as its energy. It is demonstrated that in this region, for 
temperatures in the range of r=400-600 MeV, and for jet energies above 20 GeV, the leading 
term in the series is dominant over the next-to-leading term and thus yields an estimate of the 
value of q. In these proceedings we discuss the details of the numerical calculation. 



1. Introduction 

In several formulations of jet modification Q], transport coefficients play a crucial role. 
These, in principle, non perturbative quantities encode the effect of the medium on the modifi- 
cation encountered by the hard jet. In the Higher- Twist (HT) approach [2|, transport coefficients 
represent the expectation values of operator products that have been factorized from the hard 
process of jet scattering and radiative emission in the medium. In the standard HT formalism, 
these are considered entirely non-perturbative and no attempt is made to evaluate them. In these 
proceedings, we outline the results of a reformulation of the HT scheme allowing for the evalu- 
ation of these coefficients in lattice gauge theory at finite temperature. The modification of the 
formalism has been described in Refs. Q. Here we discuss the pertinent numerical results. 

2. Formulation in Minkowski Space-Time 

Consider a single quark with a large negative light-cone momentum q~ propagating through 
a box of length L, held at a fixed temperature T = 1//3, undergoing a single scattering which im- 
parts transverse momentum to the quark. The quark then exits the medium and is observed. The 
length of the medium will be considered to be short enough such that the quark does not radiate 
during traversal. The coupling of the jet with the medium will be considered to be small enough 
such that secondary scattering is minimal. Dividing the square of the transverse momentum 
gained by L, we obtain the non-perturbative definition of the transport coefficient q, i.e., 

, 47T 2 a s f dy-d 2 y ± 2 -i*L y - + £ y> e^ E " _ . 

<?=— ; rr~ d kj_e ^ ■ ±± (n\—- F > ± (y ,y±)F ± (0)\n). (1) 

In the equation above, y~ is a light-cone separation, F 1 ™ = t a F ativ is the gauge field strength, \n) 
is a state of the thermal ensemble with a partition function Z, as is the strong coupling constant, 
and N c is the dimensionality of the fundamental representation. 
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Compared to the large energy of the hard parton, q~Q (where Q denotes the hardest scale 
in the problem), the transverse momentum imparted to the parton (k ± ) is considered to be rather 
small, k± ~ AQ, where A <k 1 is a small dimensionless constant. As a result, the (+)-component 
of the the momentum transferred, k + ~ k\/{2q~) is of the order of A 2 Q. As argued in Ref. 0, 
this operator product, q, defined over values of k + - k\/(2q~) may be expressed as the imaginary 
part of the more generalized operator product, 

. _ An 2 a s C d 4 yd 4 k j k . y 2(g-) 2 <M|F + ^(0)F+ (y)|M) 

N c J (2;r) 4 6 ' y/2q- (q + k) 2 + ie ' U) 

As defined above, Q has a branch cut in the regime where q + ~ A 2 Q, the discontinuity across the 
cut yields q. One may then define the following integral in the deep Euclidean region, 

dq + Q{q + ) 
2ni (q + + q~) 

where, the contour is taken around the point q + - -q~. We now evaluate the integral using the 
method of residues, at q + = -q~. Since q~ is very large compared to k, one may expand out the 
denominator, and the y~ dependence in the numerator to obtain, 
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The above equation represents a series of local operators, containing rising powers of derivatives 
in the numerator along with similarly rising powers of the hard scale q~ in the denominator. In 
any evaluation of I\, the presence of a large energy scale q~ will allow one to terminate the series 
at a reasonably small value of n while encumbering a minimal amount of error. 

We can now deform the contour and evaluate it over the branch cut from q + > -A 2 Q to 
q + — > oo. This yields, 
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dq I + dq + V(q + ). (5) 
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The second term in the equation above, refers to the contribution to the operator above from 
vacuum gluon radiation, i.e., the Bremsstrahlung radiation of gluons from an off-shell quark. As 
such, it contributes only in the region where the virtuality of the incoming quark is time-like 
and is independent of the temperature of the medium. Thus when T is varied, the second term 
above is a constant, while the first depends on the temperature of the medium. Comparing the 
two equations Q and PJ, subtracting the pure vacuum contributions, one may evaluate q. In the 
remainder, we discuss the the evaluation of the operators in Eq. (Hh. 



3. Euclidean space and results of calculations 

All the operator products in Eq. Q represent local operators which may be easily evaluated 
by rotation to Euclidean space. Denoting a generic Minkowski space operator product as £) > (f) = 
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Figure 1: Left panel: The temperature dependence of the local operator {F + 'F + ') scaled by 7" 4 to make it dimensionless. 
The lattice spacing is set using a non-perturbative RG prescription. Right panel: Unsealed expectation of lattice-size- 
independent correlator 2i=l,2 a f (F 31 F 3i - F 4I F 4I )/2 at finite temperature (red squares), versus expectation in vacuum 
(black circles), as a function of y3 = 4/g 2 (g is the bare lattice coupling). The plot is for n, = 5, n s = 20. 



2„(«|e"" /JH C)i(f)(32(0)|n}, we obtain the simple relation that, D > (t = 0) = i N 'A(r = 0). Where, 
N, denotes the number of time derivatives in D > (t). Using the above relation, local operator 
products in Minkowski space may be obtained from local operator products in Euclidean space. 

We report results on a (4 x n f ) 3 x n, lattice where n, is varied from 3 to 6. All calculations have 
been done with 5000 heat bath sweeps for each data point. We have used the Wilson gauge action 
for SU(2) JUIF). The scale is set on the lattice using two different renormalization group (RG) 
formulas: The first is based on the two loop perturbative RG equation for the string tension |4l[5], 
which yields the following formula for the lattice spacing, 

1 l\\g 2 \^ I \2n 2 \ 

fl i = 7-krr ex P -7TT > (6) 
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where, g represents the bare lattice coupling and represents the one dimension-full parameter 
on the lattice. Comparing with the vacuum string tension, we have used Al = 5.3 MeV. For a 
lattice at finite temperature or one with n, <K n s , the temperature is obtained as T = l/(n,a L ). The 
results for the field-strength-field-strength correlation (FFC), 2L=i,2(^ 3 '^ 3 ' _ F 4 'F 4 ')/2 with this 
choice of formula for the lattice spacing are presented in the left panel of Fig [2] The resulting 
correlation is scaled by T 4 , with T as obtained from the formula above. 

We have also set the scale using a non-perturbative approach as outlined in Ref. [6|, where 
the formula for the lattice spacing is expressed as the product of that obtained from Eq. (|6|l and 
a non-perturbative function A(g 2 ) which has been dialed to ensure that T c /Al is independent of 
g 2 . The results for the FFC, Z/=i,2(^ 3 '^ 3 ' _ F 4 'F 4 ')/2 with this next choice of formula for the 
lattice spacing are presented in the left panel of Fig[T] While the scaling with lattice size is seen 
to be much better with the non-perturbative RG equation, for temperatures above T = 400 MeV, 
both calculations yield the same value of the leading operator product in the series for I t . In 
the following we will only focus on the results above a T — 400 MeV. As a result, we will use 
the perturbative RG equation to set the scale in the evaluation of the next-to-leading terms. The 
field-strength-field-strength correlation, plotted in the left panels of Fig. [T]and Fig. [2] contain 

3 



their vacuum parts, these have not been subtracted out. In the right panel of Fig.[T[ we plot the 
unsealed correlator as a function of the coupling on the lattice and also plot the vacuum piece. 
For the points above the transition, the vacuum contribution is seen to be a small correction, and 
thus was ignored in all other plots. 
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Figure 2: Left panel: The temperature dependence of the local operator (F + 'F +I ) scaled by 7"* to make it dimensionless. 
Right panel: temperature dependence of absolute values of the local operator (F + 'F +I ) and the next-to-leading operator 
products <[F + 'i£> 4 F +, ']> [red diamond (dashed)] and ([F +i i£>*F +i ] /cf) [ green diamond (dot dashed)] . 

The next-to-leading operator product in the series for is plotted in the right panel of Fig. [2] 
Note, this operator has an extra temporal derivative and suppressed by a factor of q . The results 
of the lattice calculation of this operator are plotted, both with and without this large factor in the 
denominator (cf is chosen to be 20 GeV). For temperatures below T - 600 MeV, this next-to- 
leading correction is found to be less than 20% of the leading term. As a result, in the range of 
temperatures from 400-600 MeV, the first term in the series expansion of /; can be used to obtain 
an estimate for q. The phenomenology of these estimates is presented in Ref. [3]; for the case 
of a SU(2) quark traversing a SU(2) plasma, we obtain a virtuality averag ed§ = 0.186GeV 2 /fm, 
at a T — 400 MeV. Naively scaling this result with the number of colors and flavors for the 
case of an SU(3) gluon traversing an SU(3) plasma with two light flavors of quarks, we obtain 
q = 1.3 - 3.3 GeV 2 /fm at a T = 400 MeV. 
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